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ABSTRACT 

The broad emission lines (BELs) of quasars and active galactic nuclei (AGNs) are important 
diagnostics of the relative abundances and overall metallicity in the gas. Here we present new 
theoretical predictions for several UV BELs. We focus specifically on the relative nitrogen 
abundance as a metallicity indicator, based on the expected secondary enrichment of nitrogen 
at metallicities Z > 0.2 Z . Among the lines we consider, N III] A1750/O III] A1664, N V 
A1240/(C IV A1549 + O VI A1034) and N V/He II A1640 are the most robust diagnostics. We 
argue, in particular, that the average N V BEL is not dominated by scattered Lya photons from 
a broad absorption line wind. We then compare our calculated line ratios with observations from 
the literature. The results support earlier claims that the gas-phase metallicities near quasars 
are typically near or several times above the solar value. We conclude that quasar activity is 
preceded by, or coeval with, an episode of rapid and extensive star formation in the surrounding 
galactic (or proto-galactic) nuclei. Chemical evolution models of these environments suggest 
that, to reach Z > Z Q in well-mixed interstellar gas, the star formation must have begun >10 8 yr 
before the observed quasar activity. 

Subject headings: line formation — galaxies: formation — galaxies: nuclei — quasars: emission lines — 
quasars: general 



1. Introduction 

Measurements of the elemental abundances 
near quasars or, more generally, active galactic nu- 
clei (AGNs) are important for several reasons. For 
example, the metal abundances affect the opaci- 
ties, kinematics and overall physical structure of 
the various emitting and absorbing regions (e.g., 
Komossa & Mathur 2001, Ferland et al. 1996, 
hereafter F96 ) . On much larger scales, metal-rich 
outflows from quasars or their host galaxies might 
be important sources of metal enrichment to the 
inter-galactic gas at early cosmological times (Fri- 
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aca & Terlevich 1998, Collin & Zahn 1999a, Sadat, 
Guiderdoni & Silk 2001). Quantifying these effects 
requires explicit abundance measurements, ideally 
at a range of distances from the active nucleus. 

Quasar abundance data also provide informa- 
tion on the chemical history of the gas and there- 
fore, indirectly, on the epoch and extent of nearby 
star formation (see Hamann & Ferland 1999, here- 
after HF99, for a general review). The metals near 
high-redshift quasars might come from some of the 
first stars forming in galactic or proto-galactic nu- 
clei after the Big Bang. The constraints on star 
formation and chemical evolution near quasars are 
therefore a valuable complement to other high red- 
shift studies that use different diagnostics or probe 
larger galactic structures (e.g., Lu et al. 1996, Pet- 
tini et al. 1999, Steidel et al. 1999, Prochaska, 
Gawiser & Wolfe 2001 and references therein). 

In this paper, we develop further the theoret- 
ical basis for deriving abundances from the rest- 
frame UV broad emission lines (BELs). We note 



that although these features form within ~1 pc of 
the central AGNs (Kaspi et al. 2000), the derived 
abundances are diagnostic of the chemical evolu- 
tion in whatever (larger) region processed the gas. 

Previous BEL studies have shown that we can- 
not simply compare a strong metal line 5 , such as 
C IV A1549, to Lya to derive the overall metal- 
licity, Z (e.g., Hamann & Ferland 1993a, HF99). 
In numerical simulations of photoionized BEL re- 
gions (BELRs), the C IV /Lya flux ratio changes 
by factors of just a few (and in a non-monotonic 
way) as the metallicity (C/H) increases by orders 
of magnitude. This insensitivity to Z is caused 
mainly by the important role of the collisionally- 
excited C IV line in cooling the gas; increas- 
ing C/H ratios, which might otherwise increase 
C IV/Lya, lead instead to lower gas temperatures 
and a nearly constant C IV/Lya flux. Another 
major drawback is that C IV and Lya can form 
in spatially distinct regions (Korista ct al. 1997a) 
and, in any given region, their relative fluxes de- 
pend sharply on both the degree of ionization and 
the shape of the ionizing spectrum (see also David- 
son & Netzer 1979, Ferland 1999). 

The best BEL abundance diagnostics involve 
ions that are not important in the cooling and 
lines that have similar excitation/emission require- 
ments (forming as much as possible in the same 
gas). 

1.1. Nitrogen and Galactic Chemical Evo- 
lution 

Nitrogen lines are especially valuable because 
the N/O abundance is known to scale with O/H 
(i.e., Z) in galactic H II regions (Shields 1976, 
Pagel & Edmunds 1981, Villa-Costas & Edmunds 
1993, van Zee, Salzer & Haynes 1998, Izotov & 
Thuan 1999). This abundance behavior is at- 
tributed to "secondary" N production, whereby 
N is synthesized from existing C and O via CNO 
burning in stars (Tinsley 1980). The net result is 
N/H oc (O/H) 2 oc Z 2 and N/O oc O/H oc Z when- 
ever secondary dominates primary N production. 
The H II region data show that these scaling re- 



5 Throughout this paper we use the notation C IV and C 
III], for example, to identify the lines C IV A1549 and C 
III] A1909. These designations are unique for the lines of 
interest here (summarized in Table 1 below). Note that 
the wavelengths listed in Table 1 are often averages over 
several transitions in the same multiplet. 



lations hold approximately for Z > 0.2 Z Q . The 
relative nitrogen abundance, e.g. N/O, can there- 
fore be diagnostic of the overall metallicity even 
when direct measures of Z (such as O/H) are not 
available. 

Galactic chemical evolution models predict that 
the N production should be delayed relative to O, 
and to a lesser extent C, in environments with high 
star formation rates and short overall enrichment 
times, <2 Gyr (HF99, Hamann & Ferland 1993b, 
hereafter HF93b). The delay is significant if the 
lifetimes of the N-producing stars are comparable 
to the time scale for major increases in Z (e.g., in 
O/H, see also Henry, Edmunds & Koppcn 2000). 
The net result is that N/H and N/O can be 
below their solar values at Z = Z Q . This offset 
is only weakly present in the H II region data 
(van Zee et al. 1998) because those data represent 
more slowly evolving environments. However, the 
offset could be substantial in quasars if their host 
environments are galactic nuclei, or dense proto- 
galactic condensations, that are undergoing rapid 
evolution (HF99). 

The gas-phase N abundance is thus given gen- 
erally by 6 , 
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where the square brackets indicate logarithmic ra- 
tios relative to solar, [a/b] = log(a/6) — log(a/6) Q , 
and q is the logarithmic offset. Values of q range 
from ~0. 0-0.1 for "slow" chemical evolution (as 



6 We prefer to use oxygen (or some equivalent "ct-capturc" 
clement like Ne, Mg or Si) in these scaling relationships 
because it is a good tracer of the overall metallicity driven 
by massive-star supernovae. The behavior of N/C, in 
particular, might be more complex, depending on the shape 
of the initial mass function and the chemical evolution time 
scale. The complication is that C/O can also increase 
(by factors up to 2-3) as the system evolves, because of 
the late-time release of C from intermediate-mass stars 
(Wheeler, Sneden & Truran 1989, HF93b). Secondary N 
production will respond to this delayed C enrichment, but 
the departures of N/H or N/O from Equations 1 & 2 should 
be small because the total C+O abundance is in all cases 
dominated by O. (We assume implicitly that secondary N 
comes from both C and O, and that the primary yields of 
these elements are not dependent on Z, c.f. Henry et al. 
2000.) 
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reflected in the H II region data) to ^0.2-0.5 dur- 
ing the rapid evolution that may occur in mas- 
sive galactic nuclei (see the enrichment model in 
HF99). 

1.2. Quasar Abundance Diagnostics 

Shields (1976) proposed that several ratios of 
UV intercombination (semi-forbidden) BELs can 
provide accurate CNO abundance ratios. The 
lines include, for example, C III] A1909, N III] 
A1750 and O III] A1664. Early results gave 
evidence for an overabundance of N relative to 
C and O, indicating Z > Z Q by the argument 
in §1.1 (see also Davidson 1977, Osmcr 1980, 
Gaskell, Wampler & Shields 1981, Uomoto 1984). 
However, those studies did not consider collisional 
deexcitation of the upper energy levels (Baldwin 
& Netzer 1978). The critical densities (where the 
rate of collisional deexcitation matches the rate of 
radiative decays) are as low as ~3 x 10 9 cm -3 for 
C III] (see Table 1 below). The observed strength 
of C III] relative to C IV led to the general belief 
that BELR densities are <3 x 10 9 cm~ 3 (Davidson 
& Netzer 1979). However, more recent line 
variability studies (Peterson 1993, Goad, O'Brien 
& Gondhalekar 1993, Ferland et al. 1992, Reichert 
ct al. 1994, Korista ct al. 1995, Peterson & 
Wandel 1999, Kaspi & Netzer 1999) indicate that 
the BELR is spatially stratified with, quite likely, 
a wide range of densities. Different lines form 
in different regions depending on their ionization 
and/or density requirements (see also Baldwin 
ct al. 1995, Korista & Goad 2000). Ferland 
et al. (1992) used variability results for NGC 
5548 to estimate a density of ^lO 11 cm -3 in the 
C IV region, well above the previous upper limit 
inferred from C III]/C IV. 

Partly to overcome this uncertainty, but also 
to invoke stronger and thus easier-to-measure 
BELs, Hamann & Ferland (1992), HF93b and F96 
developed N V A1240/C IV and N V/He II A1640 
as abundance diagnostics. Their results indicate 
that enhanced N abundances and solar or higher 
metallicities are typical of BELRs (see also Osmer 
ct al. 1994, Laor ct al. 1994 and 1995, Korista et 
al. 1998, Dietrich et al. 1999, Dietrich & Wilhelm- 
Erkens 2000, Constantin et al. 2001, Warner et al. 
2001). 

Independent confirmation of the BEL abun- 
dances has come from quasar associated absorp- 



tion lines (AALs), which also reveal typically high 
metallicities (Z > Z Q ) and probably enhanced 
N/C ratios in gas near the quasars (Wampler et al. 
1993, Petitjean, Rauch & Carswell 1994, Hamann 
1997, Petitjean & Srianand 1999, HF99). The 
overall implication is that quasar environments arc 
metal-rich and thus chemically "mature" — hav- 
ing already undergone substantial star formation 
•prior to the observed quasar epochs. This result 
appears to hold generally, even for quasars at the 
highest measured redshifts (HF99, Constantin et 
al. 2001, Warner et al. 2001). 

Here we test and further quantify the abun- 
dance sensitivity of various BEL ratios. Through- 
out this paper, we define solar abundances by 
the meteoritic results in Grevessse & Anders 
(1989). We employ atomic data from Verner, 
Verner & Ferland (1996) and from the compila- 
tion given in HAZY — the documentation sup- 
porting the spectral synthesis code, CLOUDY, 
which we also use below for line flux calculations 
(Ferland et al. 1998). CLOUDY and HAZY 
are both freely available on the World Wide Web 
(http:/ /www.pa.uky.edu/^gary/cloudy). 

2. Analytic Estimates of the Line Flux 
Ratios 

Analytic estimates of the line flux ratios can 
yield both physical insight and accurate abun- 
dance results. The theoretical ratio of two col- 
lisionally excited lines emitted from the same vol- 
ume by two-level atoms is, 

Fx n n u 1^2312 / A£2 - A£i \ 

— = Q — ^— : exp — (3 

F 2 ni2 T 2 Ai.g n \ kT e J 

where Q is defined by 

1 J Ik 

Q - Y~^t ( 4 ) 

' n cr i 0! 

Stimulated emission and interactions with the con- 
tinuum radiation field are assumed to be negligi- 
ble. For each line designated 1 or 2, F is the mea- 
sured flux, A is the rest wavelength, AE is the 
photon energy, T is the energy-averaged collision 
strength, n; and gi are the number density and sta- 
tistical weight of the lower energy state, (3 is the 
photon escape probability (0 < (3 < 1), and 
the critical density of the upper energy state. T e 
is the electron temperature and n e is the electron 
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density. Table 1 lists critical densities for the lines 
of interest here. 

The factor Q in Equation 3 accounts for possi- 
ble photon trapping and collisional deexcitation. 
If the densities are low, such that n e -C n cr (3 
for both lines in a given ratio, then Q « 1 and 
collisional deexcitation is not important. If the 
densities are high such that n e 3> n cr j3 for both 
lines, then collisional deexcitation is important 
and Q « n cr ij3i/n cr 2(32- (The level populations 
approach LTE in this limit.) If the line photons 
escape freely (e.g., for forbidden or semi- forbidden 
lines with low oscillator strengths) , then /?i , (3 2 « 1 
and the correction factor at high densities is sim- 
ply Q « n cr i/n cr2 . 

The lines become "thermalized" and thus lose 
their abundance sensitivities 7 only in the limit 
where j3 <C 1 and n e » n cr (3, that is, where 
the lines are optically thick and the density ex- 
ceeds the "effective" critical density of n cr /3. The 
abundance sensitivity remains in all other circum- 
stances, namely, (1) in all cases where the lines are 
optically thin - including the high density limit 
(where n e ^> n cr (3), and (2) at large line opti- 
cal depths if the density is low enough to satisfy 
n e <C n cr (3 (see also HF99). 

If the ions are primarily in their ground states 
(or ground multiplets, such that ni approximates 
the total ionic density), then we can rewrite 
Equation 3 as, 

Fi nn f(X\) fXA fC\ 

where T4 is the gas temperature in units of 10 4 K, 
f{X\) is the fraction of clement X\ in ion stage 
X{, etc., and X1/X2 is the abundance ratio 
by number. B and C contain the various rate 
coefficients and physical constants from Equation 
3. Table 2 provides values of B and C for the line 
ratios of interest here. B has a weak temperature 
dependence (leftover from the ratio of collision 
strengths) that is negligible for our applications 
(c.f.. Netzer 1997). The tabulated B values 

7 We can see that the abundance sensitivity is lost in 
the high density, optically thick limit because Q — > 
n-cri/3i/"-cr2/32 tx Y2T2/Y1T1 oc T 2 n( 2 /Tin(i, where r is 
the line optical depth and (3 ~ 1/r for r ^> 1 (Frisch 1984). 
Plugging into Equation 3 shows that the line ratios do not 
depend on the relative abundance in this regime. 



assume T e — 10 4 K. Table 2 also lists the flux 
ratios derived from Equation 5 in the two limiting 
optically thin cases, Q — 1 and Q = n cr i/n cr 2, 
assuming f{X\)/ f{X 3 2 ) = 1 and solar abundances. 
The temperatures used for these flux ratios are 
based on the numerical simulations in §3.1 below, 
namely, T 4 = 0.8 for N II] /C II], T 4 = 2.2 for the 
ratios involving N V, and T4 = 1.5 for all others. 
Note that the Q = n cr i/n cr 2 flux ratios listed in 
Table 2 for the permitted lines are not expected 
to occur in BELRs, because these lines are very 
optically thick (§3.3.3, Figure 6 in Hamann et al. 
1995). 

3. Numerical Simulations 

Modern numerical simulations yield more 
reliable abundance results because they self- 
consistently treat the ionization, temperature, 
radiative transfer and line emission in realis- 
tic BELR clouds. We use the numerical code 
CLOUDY, version 96.00 (Ferland et al. 1998) to 
predict BEL fluxes in different circumstances. We 
assume the emitting clouds are free of dust and 
photoionized by the quasar continuum radiation. 
The degree of ionization depends mainly on the 
shape of the incident spectrum and the ionization 
parameter, U = ^n/cnn, where «h is the vol- 
ume density of hydrogen particles (H I+H II) and 
<E>h (units = cm~ 2 s _1 ) is the flux of hydrogen- 
ionizing photons incident on the clouds. The to- 
tal column density in all of our calculations is 
log7V H (cm- 2 ) = 23.5 (c.f. Ferland & Persson 
1989). 

3.1. A Typical Cloud Structure 

Figure 1 illustrates the temperature and ion- 
ization structure of a "typical" BELR cloud hav- 
ing U = 0.1, riH = 10 10 cm~ 3 , solar abundances, 
and an incident spectrum defined by Mathews & 
Ferland (1987, hereafter MF87). We consider the 
effects of other continuum shapes and cloud prop- 
erties below. The model cloud in Figure 1 pro- 
duces strong emission in most of the lines of in- 
terest here (see §3.2 below). The main exceptions 
are C II] and C III], for which n H = 10 10 cm~ 3 
is above their critical densities (Tabic 1). Pairs of 
ions that overlap substantially in Figure 1 are good 
candidates for abundance diagnostics. The figure 
shows, for example, that N+ 4 resides within the 
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He ++ zone, that the C +3 and +5 zones bracket 
the N+ 4 region, and that C +2 -N+ 2 and C +3 -N+ 3 
define cospatial zones that both overlap with +2 
(see also Netzer 1997). 

Figure 2 shows more specifically where the dif- 
ferent lines form within our "typical" BELR cloud 
(as defined above for Figure 1). In particular, the 
plot shows the line emissivities, J, multiplied by 
the local line escape probabilities, (3, as a function 
of spatial depth into the cloud. The emissivities 
are further multiplied by the spatial depth, D, in 
Figure 2 to offset the tendency for the logarith- 
mic scale to exaggerate the emission from small 
depths. Lines whose emission regions significantly 
overlap in Figure 2 may be good abundance diag- 
nostics (see also F96). 

3.2. Predicted Line Strengths and Ratios 

Figures 3 and 4 show theoretical BEL equiva- 
lent widths and flux ratios, respectively, emitted 
by clouds with different nn and $h- Other input 
parameters are the same as Figure 1. (See Korista 
ctal. 1997a, Baldwin et al. 1995 and Ferland 1999 
for many more plots similar to Figure 3.) 

We next consider the effects of different contin- 
uum shapes and metallicities ranging from Z = 0.2 
Zq to 10 Z Q . Table 3 lists the abundances used 
for each Z (see also the Appendix). The metals 
are all scaled keeping solar proportions, except for 
nitrogen, which is scaled according to Equation 
1 with 5 = 0. Note that by choosing q = we 
will derive maximum nitrogen line strengths for a 
given Z ', leading to conservatively low estimates of 
Z when compared to the measured line ratios (§7 
below). The change in relative helium abundance 
is small and therefore unimportant. We scale He 
from solar such that the change in He mass frac- 
tion, AY", equals the change in mctallicity, i.e., 
AY/AZ = 1 (Baldwin et al. 1991, HF93b and 
references therein). 

If BELRs include simultaneously wide ranges 
in fin and <I>h, then each line will form primar- 
ily in regions that most favor its emission. This 
situation has been dubbed the Locally Optimally- 
emitting Cloud (LOC) model (Baldwin et al. 
1995). It is consistent with observational re- 
sults and a natural extension of multi-zone mod- 
els (Rees, Netzer & Ferland 1989, Peterson 1993, 
Brotherton et al. 1994a, Baldwin et al. 1996, 



Hamann et al. 1998). It also has the tremen- 
dous advantage of not requiring specific knowl- 
edge of jih and $h- The total line emission is 
simply an integral over the iir-^h plane (e.g., 
Figure 3), with weighting factors that specify the 
relative numbers of clouds with high versus low 
nn and high versus low $h (see also Ferguson 
et al. 1997). We adopt equal weighting per 
decade in the nn-^H distribution, which is al- 
ready known to provide a good match to typi- 
cal AGN spectra (Baldwin 1997, Korista & Goad 
2000). The full parameter ranges included in the 
LOC integration are 7 < lognH(cm~ 3 ) < 14 and 
17 < log$ H (cm- 2 s _1 ) < 24. 

Figure 5 shows predicted line flux ratios as a 
function of Z for three different incident spectra in 
the LOC calculations 8 . The incident continua are 
i) MF87, as in Figures 1-4, ii) a "hard" power law 
with index a = — 1.0 (/„ oc v a ) across the infrared 
through X-rays, and iii) a segmented power law 
that approximates recent observations, namely, 
a = —0.9 from 50 keV to 1 keV , a = —1.6 from 1 
kcV to 912 A, and a = -0.6 from 912 A to 1 /im 
(Zheng et al. 1996, Laor et al. 1997). All of the 
model spectra have a steep decline at wavelengths 
longer than 1-10 microns to prevent significant 
free-free heating in dense clouds (Ferland et al. 
1992). The true spectral shape is not well known. 
The MF87 spectrum remains a good guess for 
the continuum incident on BELRs, in spite of 
its differences with recent observations (Korista, 
Ferland & Baldwin 1997b). In any event, the 
three spectra used for Figure 5 span a wide range 
of possibilities — from a strong "big blue bump" 
in the MF87 continuum to none at all in the 
a = —1.0 power law. Figure 5 shows that the 
uncertainty in the continuum shape has little effect 
on the line ratios considered here (except, perhaps, 
for N V/He II, which is discussed extensively by 
F96 and HF99, see also §3.3.4 and §4 below). 

It is important to note that the results in Figure 
5 do not depend on the validity of the LOC model. 
Most of the line ratios, notably N III]/0 III], do 
not vary sharply with nn or <E>h in the regions 
where the lines individually are strong (Figures 3 
and 4). The LOC integrations are therefore simply 

The results in Figure 5 are not sensitive to the LOC 
integration limits. For example, revising the limits to 
8 < logn H (cm- 3 ) < 12 and 18 < log ^(cm^s^ 1 ) < 24 
changes the line ratios by <0.1 dex at a given Z. 
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a convenient way of summarizing the information 
in plots like Figure 4. For line ratios that do have 
significant tih-^h sensitivities, such as N V/C IV, 
the LOC results should minimally represent the 
average properties of diverse quasar samples. 

3.3. Parameter Sensitivities 

3.3.1. Density and Ionizing Flux 

Most of the lines are weak near the lower 
right and upper left extremes in Figure 3 because 
the ionization parameter is too low or too high, 
respectively. Some lines are also weak in the upper 
right because, for any given U, their emission 
can be suppressed at high nn (§2, Korista et al. 
1997a). 

Figure 4 shows the sensitivity of the line ratios 
to riH and $h- We discuss the implications for 
abundance estimates in §4 below. 

3.3.2. Turbulence and Continuum Pumping 

One poorly known parameter in BELR clouds 
is the internal Doppler velocity vd- Our calcu- 
lations assume vd is strictly thermal (i.e., vd ~ 
v t h = y/2kT/m). Larger Doppler velocities would 
lead to lower line optical depths and reduced ther- 
malization effects at high densities (§2 and §3.3.3). 
Larger vd can also enhance the line emission via 
"continuum fluorescence," that is, resonant ab- 
sorption of continuum photons followed by ra- 
diative decays. If the absorption lines are opti- 
cally thick, then this fluorescent contribution to 
the BELs will depend much more on vd than the 
abundances (e.g., Ferguson, Fcrland & Pradhan 
1995, Hamann et al. 1998). 

We examine these effects by varying vjj from 
to 3000 km s _1 in the fiducial cloud model 
portrayed in Figure 1. Bottorff et al. (2000) 
present more extensive calculations using different 
cloud parameters. The overall result is that 
continuum pumping is insignificant in C IV and all 
of the intercombination lines for any vd < 3000 
km s _1 . (In C IV the collisional contribution is 
overwhelming; in the intercombination lines the 
optical depths arc too low at large vjj.) Pumping 
is similarly not important for the other lines 
considered here as long as vd < 300 km s" 1 . 
Pumping begins to exceed collisional excitation 9 

'Note that continuum pumping can be relatively more 



in N III, O III and O VI only for v D > 300 
to 800 km s _1 , and in N V and C III for 
vrj > 1000 km s _1 . The relevant line ratios are 
even less sensitive to Aw than the line strengths 
individually. For example, N V/(C IV+O VI) 
changes by <0.1 dex across the range < Av < 
1000 km s _1 . Since the actual value of vd in 
BELRs is unknown (Bottorff et al. 2000), we will 
assume that vd is not large enough to affect the 
line ratios of interest here. 

3.3.3. Line Optical Depths 

Our calculations indicate that the intercombi- 
nation lines are all optically thin, and the per- 
mitted lines are all optically thick, in the regions 
where their emission is strongest (see also Ferland 
et al. 1992, Hamann et al. 1998). The surprising 
lack of density dependence in the N IV] /C IV ra- 
tio in Figure 5 (c.f.. Table 2) is caused by the fact 
that (3n cr for N IV] is similar to (3n cr for C IV, 
so the lines behave similarly with density. This 
behavior depends on C IV becoming thermalized 
at high densities, where it no longer has a direct 
sensitivity to the abundances (§2). Thermaliza- 
tion can occur in other lines as well, most notably 
C III, but also in some of the intercombination 
lines. However, thermalization in limited parts of 
the riH-$H plane is not a problem for abundance 
studies if BELRs have, as expected, a range of 
physical conditions (Ferland et al. 1992, Broth- 
crton et al. 1994a, Baldwin et al. 1995, Bald- 
win et al. 1996). Each line's emission should be 
dominated by clouds where thermalization is not 
important. 

3.3.4- Metallicity and Abundance Ratios 

The predicted line flux ratios in Figure 5 scale 
almost linearly with the abundance ratios (N/O 
and N/C) and therefore also with the metallicity 
(because Z oc N/O oc N/C by Equation 2). Over- 
all, there is remarkably good agreement between 
the numerical and analytic results. For example, 
the numerical results at Z = Z Q lie generally be- 
tween the analytic ratios derived for Q = 1 and 
Q = n cr i/n cr2 in Table 2. 

important for a given vu in environments with lower riH 
and TVjj , such as the narrow emission line regions of AGNs 
(Ferguson et al. 1995, Nctzcr 1997). 
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It is interesting to note that N V/He II also 
scales roughly linearly with Z in Figure 5, even 
though N/He oc Z 2 (§1). This behavior is tied to 
the energy balance in BELR clouds; higher metal- 
licities increase the efficiency of metal-line cooling 
and thus lower the temperature. The temperature 
drops by exactly the amount needed to maintain 
the total line emission and preserve the overall en- 
ergy balance. N V/He II reacts strongly to tem- 
perature changes because it compares a collision- 
ally excited line to a recombination line (c.f. Equa- 
tions 2, 3 and 6 in HF99). The rise in N V/He II 
with increasing N/He is therefore moderated by 
the dropping temperature. The net result is that 
N V/He II scales roughly like N/O oc Z (i.e., like 
the enhancement of N relative to the other metal 
coolants; see also F96, Ferland et al. 1998 and 
HF99). 

4. Preferred Line Ratio Diagnostics 

The theoretical accuracy of any line ratio abun- 
dance diagnostic depends on a variety of factors, 
such as 1) the temperature sensitivity of the ratio 
(Equation 5), 2) the similarity of the ionization 
potentials and critical densities, 3) the extent to 
which the line-emitting regions overlap spatially 
(Figures 1-3), and 4) the constancy of the ratio 
across a range of plausible Uh and &h (Figure 4). 

Based on these criteria, the most reliable abun- 
dance diagnostic among the intercombination lines 
is N III]/0 III]. Note that +2 can occupy a larger 
region than N +2 (Figure 1), and therefore low col- 
umn density clouds might contribute significant 

III] emission but little N III]. Comparing ob- 
served values of N III]/0 III] to our calculations 
should therefore yield, in the worst case, lower lim- 
its to N/O and thus Z. 

N III]/C III] is also useful, but it may be less 
robust than N III]/0 III] because it has a stronger 
temperature dependence and the relevant critical 
densities differ by a factor of nearly 2 (see Tables 

1 and 2). The net result is a greater sensitivity 
to non-abundance parameters (e.g. Figure 4). For 
example, a larger/smaller proportion of low den- 
sity clouds in the LOC integration would have led 
to ^2 times lower/higher predicted N III]/C III] 
ratios in Figure 5. We also recall (footnote 6) that 
N/C may not be as clearly tied to the metallicity 
as N/O.' 



Among the collisionally-excited permitted lines, 
the most robust diagnostics should be N III/C III 
and N V/(C IV+O VI) because of their relatively 
weak temperature dependencies. Note, in partic- 
ular, that the ionization and excitation energies 
of C+ 3 and 0+ 5 bracket those of N+ 4 along the 
lithium isoelectronic sequence. The combined ra- 
tio N V/(C IV+O VI) therefore has lower sensi- 
tivities to temperature, ionization and continuum 
shape compared to the separate ratios N V/C IV 
and N V/O VI. 

N V/He II is sensitive to the temperature 
and ionizing continuum shape (§3.3.4) and, like 
N III]/0 III], it may under represent the N/He 
abundance because N+ 4 can occupy a smaller 
region than He ++ (inside the He ++ zone, F96, 
Figure 1). Nonetheless, HF93b and F96 showed 
that by adopting BELR parameters that maximize 
the predicted N V/He II (for example the a = 
— 1.0 power law in Figure 5), this ratio can yield 
firm lower limits on N/He and therefore Z. 

5. Is N V Enhanced by Lya? 

Several studies have noted that the N V emis- 
sion line might be enhanced by Lya photons scat- 
tered in a broad absorption line (BAL) wind (Sur- 
dcj & Hutsemekers 1987, Turnshek 1988, Hamann, 
Korista & Morris 1993, Turnshek et al. 1996, Kro- 
lik & Voit 1998). The underlying premise is that 
all quasars have BAL outflows, but they are de- 
tected only when the gas happens to lie along our 
line of sight to the continuum source. The BAL 
gas is believed to be located outside the radius 
of both the continuum emission region and the 
BELR (see, for example, Figure 4 in Cohen et 
al. 1995). If BAL winds include N+ 4 ions out- 
flowing at ^5900 km s _1 (corresponding to the 
wavelength difference between Lya and N V), then 
these ions can resonantly scatter Lya BEL (and 
continuum) photons into our line of sight - thus in- 
creasing the observed N V emission strength. This 
scattering contribution from a BAL wind is not in- 
cluded in our calculations (§2 and §3). 

However, there is strong evidence that the scat- 
tering contributions to N V BELs are typically 
small. In particular, the average N V BAL inter- 
cepts only ~30% of the incident Lya flux (based 
the depth of the average N V BAL trough ex- 
trapolated across Lya; Hamann & Korista 1996). 
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Therefore, only ^30% of the Lya photons incident 
on BAL regions are available for scattering into the 
N V BEL. (Krolik & Voit's [1998] estimate of the 
scattered N V flux is at least ~3 times too large 
because they assumed that N V BALs intercept 
all of the incident Lya photons.) 

The actual scattering contributions to measured 
BELs depend further on 1) the scattered line pro- 
files (i.e., how the scattered line flux is distributed 
in velocity), and 2) the global covering factor of 
the BAL wind (i.e., the fraction of the sky cov- 
ered by BAL gas as seen by the central continuum 
source). The total scattered flux scales linearly 
with the covering factor. Detailed line scatter- 
ing calculations assuming an average BAL velocity 
profile and a nominal 12% covering factor yield 
total scattered N V fluxes that are only ^25% 
of the average N V BEL measured in quasars 
(Hamann & Korista 1996). The 12% covering fac- 
tor is based on the ~12% detection frequency of 
BALs in quasar samples (Foltz et al. 1989 and 
1990). There is some evidence from polarization 
studies that the average covering factor might be 
20% to 30% (Schmidt & Hines 1999). However, 
30% is probably a firm upper limit because larger 
covering factors would lead to too much scattered 
flux filling in the bottoms of observed BAL troughs 
(Hamann et al. 1993). 

The calculations by Hamann & Korista (1996) 
also show that the scattered line flux from a BAL 
wind is likely to appear over a wider range of ve- 
locities than typical BELs. For scattering to domi- 
nate the measured N V BEL, the N+ 4 ions in BAL 
winds, with velocities reaching typically 10,000- 
20,000 km s , must scatter Lya plus continuum 
photons into BEL profiles that have average half- 
widths of only ^2500 km s _1 (e.g., Brotherton 
et al. 1994b). Isotropic scattering of the Lya 
photons can therefore be ruled out, because that 
would produce flat-topped N V emission profiles 
with half- widths of ~5900 km s _1 . Scattered con- 
tinuum photons would appear over an even wider 
range of line velocities (±20,000 km s _1 based on 
the average N V and C IV BAL troughs, Hamann 
& Korista 1996). Highly anisotropic wind geome- 
tries are required, but that causes other problems 
if the scattered flux dominates, namely, BELs that 
should be broader in quasars with observed BALs 
compared to quasars without BALs (see Hamann 
& Korista 1996 and Hamann et al. 1993 for spe- 



cific calculations). These differences are not ob- 
served in quasar BELs (Weymann et al. 1991). 

Another important point is that observed N V 
emission profiles closely resemble those of other 
BELs (e.g., C IV) measured in the same spectrum 
(Osmeretal. 1994, Laor et al. 1995, Constantin et 
al. 2001). This occurs naturally for lines forming 
in more or less the same BELR gas. However, the 
profile similarities become a serious problem if N V 
is attributed to scattering in a high- velocity wind, 
while the other BELs are collisionally excited in a 
separate region (the usual BELR) whose velocity 
field is believed to be smaller and mostly not radial 
(Tiirlcr & Courvoisier 1997, Korista et al. 1995, 
Peterson & Wandel 1999). 

Finally, reverberation studies indicate that high 
ionization BELs such as N V and He II form 
generally closer to the continuum source than Lya 
and other lower ionization lines (Korista et al. 
1995, Dietrich & Kollatschny 1995, Wanders et 
al. 1997, O'Brien et al. 1998). In contrast, N V 
BALs are observed to suppress the Lya emission 
flux, presumably because the N V BAL region 
lies outside the Lya BELR radius. These results 
are apparently contradictory if N V BELs are 
attributed to scattering in a BAL wind. 

All of these difficulties suggest that measured 
N V BELs are controlled by normal BELR emis- 
sion, and not by scattered flux from a BAL out- 
flow. 

6. Observed Line Ratios 

Figure 5 includes observed BEL flux ratios from 
several studies (see figure caption) . The data rep- 
resent quasars over a range of redshifts and lu- 
minosities. For example, Dietrich & Wilhelm- 
Erkens (2000) measured luminous quasars at red- 
shifts z > 3, while Laor et al. (1994, 1995) and 
Wills et al. (1995) observed relatively less lumi- 
nous sources at z < 1. The data plotted from 
Boyle (1990) derive from a mean spectrum that 
emphasizes relatively low luminosity quasars with 
z ~ 1-2. 

We include data points from the literature only 
if both lines in a given ratio were measured. This 
criterion selects against weak lines and might bias 
the plotted data toward, for example, smaller 
N V/He II ratios or larger N III]/C III]. However, 
inspection of the original references indicates that 
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any such biasing is not severe. In particular, if 
we include upper limits on undetected weak lines, 
the resulting ratios are still consistent with the 
ranges shown for directly measured values. Ratios 
involving two weak lines, such as N III]/0 III], 
should also be unbiased because such lines are 
generally either both present or both absent in 
the data. In addition, the measurements from the 
mean spectrum in Boyle (1990) are not subject 
to these selection biases, yet they are similar to 
the other average values in Figure 5. We therefore 
expect the data in Figure 5 to be representative of 
the original quasar samples. 

Several of the observed BELs can be blended 
with adjacent features. We use "deblended" 
measurements in all cases except for C III], which 
can be blended with Si III] A 1892 and Al III] 
A1859. Efforts to deblend the C III] complex 
(Laor et al. 1995, Steidel & Sargent 1991) have 
shown that C III] contributes typically 65% to 
90% to the total flux. If the original study made 
no attempt to deblend these lines, we multiplied 
the flux reported for the overall blend by 0.8 to 
approximate the C III] contribution alone. 

The observed value of N IV] /O IV] in Figure 5 
is based on a single N IV] /O IV] measurement by 
Laorctal. (1995). This measurement is unreliable 
because N IV] is weakly detected and O IV] is 
severely blended with Si IV AA1394,1403 emission. 

7. Abundance Results 

The observed line ratios are shown in Figure 
5 on the left-hand side of each panel for conve- 
nience. The metallicities implied these data are 
determined by noting the value of Z where the 
observed ratios intersect the theoretical curves. 
We conclude from these comparisons that quasar 
metallicities are typically near or a few times 
above the solar value. In particular, the most 
robust diagnostics, such as N III]/0 III] and 
N V/(C IV+O VI) (§4), indicate typically Z ~ 1- 
3 Z Q . C III]/N III] and C III/N III provide similar 
results. N V/He II also strongly supports Z > Z Q , 
and is consistent with Z ~ 1-3 Z Q if we adopt the 
relatively "hard" power law continuum, a = — 1. 
N IV] /O III] is less reliable theoretically (§2 and 
§3.1) and shows considerable scatter among the 
observations (see also Warner et al. 2001). We 
therefore do not consider this ratio further. We 



also disregard the other N IV] line ratios because 
they have even larger uncertainties (§4 and §6). 

Actual quasar metallicities could be as much 
as 2-3 times higher than these estimates because 
the theoretical calculations (Table 2 and Figures 
4 and 5) are based on q = in Equation 1. In 
fact, q ~ 0.2-0.5 may be more appropriate for the 
rapid chemical enrichment occurring near quasars 
in galactic nuclei (HF93b, §1.1). The results 
presented here are therefore broadly consistent 
with quasar metallicities 10 being typically in the 
range Z <~ 1-9 Z Q . 

8. Discussion 

It is important to keep in mind that the metal- 
licities discussed here are based on the assump- 
tions of secondary nitrogen production and N oc 
Z 2 scaling (§1.1). Actual abundance ratios might 
vary if, for example, there was a recent short- 
duration starburst (which could magnify differ- 
ences in the time-dependent yields of C, N and 
O; Lu, Sargent & Barlow 1996, Coziol et al. 1999, 
Henry et al. 2000). It is also possible that the 
scaling relations in Equations 1 and 2 break down 
at Z > Z if the yields of C and O depend signifi- 
cantly on Z (Henry et al. 2000). In addition, rare 
quasars with very strong nitrogen lines (Osmer 
1980, Baldwin et al. 2001) might not have metal- 
licities as high as Figure 5 would suggests; their 
environments might be unusually nitrogen rich be- 
cause of localized "pollution" by Wolf-Rayet stars, 
luminous blue variables or rare planetary nebulae 
(Davidson et al. 1986, Kobulnicky et al. 1997, 
Tajitsu et al. 1999). 

Nonetheless, an approximate N oc Z 2 scal- 
ing seems well established for most circumstances 
(§1.1), and so the metallicities inferred from Fig- 
ure 5 should be approximately correct. It is also 
reassuring that the results derived here from the 
BELs are in good agreement with mctallicity esti- 
mates based on intrinsic absorption lines (§1.2), 
which do not in any way rely on assumptions 
about the nitrogen production. The metallicities 
implied by all these diagnostics are typically near 
or a few times above the solar value. 



Some previous BEL studies favored the high end of this 
metallicity range (e.g., HF93b, Fcrland et al. 1996, Korista 
et al. 1998, Dietrich et al. 1999) largely because they 
assumed q > for the nitrogen enrichment. 
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If these gas-phase metallicities result from the 
normal evolution of surrounding stellar popula- 
tions, then those populations must already be 
fairly evolved at the observed quasar epoch (but 
see Collin & Zahn 1999b for an alternative hypoth- 
esis). In particular, chemical evolution models us- 
ing "standard" initial mass functions show that 
Z > Zq occurs in well-mixed interstellar gas only 
after most of the original gas has been converted 
into stars and stellar remnants (e.g., HF93b). At 
high rcdshifts, this evolution must have occurred 
quickly (with a high star formation rate) to achieve 
Z > Z Q by the time the quasars became observ- 
able. For example, at redshift z ~ 5 the available 
evolution time is only ~l-2 Gyr, depending on 
cosmological models (see Figure 1 in HF99). 

Such rapid evolution may seem extreme by solar 
neighborhood standards, but it is well within the 
parameters derived independently for dense proto- 
galactic condensations (Ostriker & Gnedin 1997, 
Cen & Ostriker 1999) and (at least the cores of) gi- 
ant elliptical galaxies (HF93b, Friaca & Terlevich 
1998, Granato et al. 2001, HF99 and references 
therein). It is also consistent with other studies 
suggesting that the degree of chemical enrichment 
in any cosmic structure depends more on its den- 
sity (i.e. the depth of its gravitational potential) 
than on its redshift (Pettini 2001). Quasars pro- 
vide direct observational evidence for high metal- 
licities and rapid evolution in dense galactic or 
proto-galactic nuclei at high redshifts (see also Di- 
etrich & Hamann 2001). Complementary evidence 
has come from observations of quasar host galax- 
ies (Nolan et al. 2001), which suggest that the 
hosts are typically massive ellipticals that have al- 
ready formed by the epoch of peak quasar activity 
[z ~ 2.5). 

Chemical enrichment models designed specifi- 
cally for dense galactic nuclei (HF93b, Friaca & 
Terlevich 1998) suggest that the time needed to 
reach Z ~ Z Q in the gas is only of order 10 8 yr. 
We conclude, therefore, that quasar host galax- 
ies begin an episode of major star formation at 
least ~10 8 yr before the onset of visible quasar 
activity. These star forming episodes might be 
among the first to occur in collapsed structures af- 
ter the Big Bang. It is clearly desirable now to ex- 
tend the quasar abundance analysis to the highest 
possible redshifts, and compare the results across 
a wide range of rcdshifts, luminosities and other 



quasar/host galaxy properties. 

9. Summary 

We present new theoretical calculations show- 
ing the abundance sensitivities of various UV 
broad emission lines in quasar/ AGN spectra. We 
focus on ratios involving C, N and O, and rely 
on the secondary enrichment of nitrogen to pro- 
vide a crude metallicity indicator. The most ro- 
bust abundance probes among the line ratios we 
consider are N III]/0 III], N V/(C IV+O VI) 
and N V/He II. Comparing our calculations to 
observed data from the literature suggests that 
quasar metallicities are typically near or several 
times higher than the solar value. These metal- 
licities probably result from the rapid evolution of 
stellar populations in galactic (or proto-galactic) 
nuclei at high rcdshifts [z > 2-3). Chemical evo- 
lution models of these environments suggest that, 
to reach Z > Z Q in the gas, the episode of major 
star formation began >10 8 yr before the quasars 
become observable. 

We are grateful to Matthias Dietrich, Hagai 
Netzer, Bassem Sabra and an anonymous referee 
for comments on this manuscript. FH acknowl- 
edges financial support from the NSF via grant 
AST 99-84040. GF thanks the NSF for sup- 
port through AST-0071180 and NASA for support 
through its LTSA program (NAG5-8212). 

APPENDIX 

Throughout this paper, we adopt the common 
usage definition of metallicity, 

Z mctals/H O/H 

Zq" ~ (metals/H) ~ (O/H) ^ ' 

where O, H and "metals" represent abundances by 
number. However, Z is strictly speaking a mass 
fraction, such that X + Y + Z = 1, where X and 
Y represent the fractions in hydrogen and helium. 
If we define 

= mctals/H 
? ~ (metals/H) 1 ' 

to simplify the notation, and adopt AY/ AZ = 1 
(§3.2), then it is easy to show that the correct 
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relationship between the mass fraction, Z, and the 
scale factor, £, is 



Z X & + Y + (2£ - 1)Z 

where Xq , Y© and Z are the Solar mass fractions 
(Xq + Yq + Zq = 1). This expression yields 
a limiting value of Z/Z Q -» 25 (Z ->■ 0.5) 
for £ — > oo. It is actually £ that we call the 
"metallicity" and refer to as "Z" throughout this 
paper (everywhere except this Appendix, e.g. in 
Table 3). For £ < 3, the approximation Z/Z & £ 
is accurate to within 10%. However, at the highest 
value of £ = 10 in Table 3, the corresponding mass 
fraction Z is just 7.4 Z©. 
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FIGURE CAPTIONS 

Figure 1. — Ionization fractions and gas tem- 
perature, T4, versus spatial depth into a "typical" 
BELR cloud (U = 0.1, n H = 10 10 cm" 3 , solar 
abundances and an incident spectrum defined by 
MF87). The temperature is indicated in the up- 
per panel by the dotted curve and by the vertical 
scale at the right. Ionization fractions are plotted 
in the upper panel for helium ions (solid curves) 
and hydrogen (dashed), and in the lower panel for 
nitrogen (solid curves), carbon (dashed) and oxy- 
gen (dash-dot). See §3.1 

Figure 2. — Depth-weighted line emissivities, 
multiplied by the escape probability (D/3J), are 
plotted against spatial depth into the same "typi- 
cal" BELR cloud described in Figure 1. The lines 
can be identified be reference to Table 1. The 
solid curves represent lines of nitrogen, the dashed 
curves represent carbon lines, the dash-dot curves 
represent oxygen lines, and the dotted curve rep- 
resents Hell. 

Figure 3. — Theoretical rest-frame equivalent 
widths (all measured relative to the continuum at 
1216 A) as a function nn and $h m BELR clouds. 
All other parameters are the same as Figure 1. 
Contours appear every 0.25 dex, with the lowest 
(bold) contour corresponding to dex and other 
solid contours marking the 1 or 2 dex boundaries. 
The ionization parameter, U, is constant along 
diagonals from the lower left to the upper right 
in each panel. 

Figure 4. — Theoretical BEL flux ratios for 
clouds with different 7Jh and $h, based on the 
same calculations as Figure 3. The bold outer 
contours define the boundary within which both 
lines in the ratio have rest-frame equivalent widths 
>1 A. Other contours are defined as in Figure 3. 

Figure 5. - Theoretical line flux ratios de- 
rived from LOC integrations are shown as a func- 
tion of the metallicity, Z/Zq, for three different 
incident continuum shapes. The solid curve in 
each panel corresponds to the MF87 spectrum, 
while the dashed and dotted curves represent the 
a = — 1.0 and segmented power laws, respectively 
(see §3.2). The symbols in each panel mark the 
average measured line ratios from different litera- 
ture sources: * = Wills et al. (1995), O = Laor 



ct al. (1994 and 1995), □ = Dietrich & Wilhelm- 
Erkens (2000), x = Boyle (1990), and A = Uo- 
moto (1984). The symbols are shown at arbitrary 
horizontal locations in each plot. The vertical bars 
associated with each symbol show the range of ra- 
tios in the middle ~2/3 of the measured distribu- 
tion (after excluding the highest ~l/6 and lowest 
~l/6 of measured values). 
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Tabic 1: Critical Densities 



Transition 


n cr (cm 3 ) 


C II] A2326 


3.16e09 


C III A977 


1.59el6 


C III] A1909 


1.03el0 


C IV A1549 


2.06el5 


N II] A2142 


9 57e09 


N III A991 


8.09el5 


N III] A1750 


1.92el0 


N IV] A1486 


5.07el0 


N V A1240 


3.47el5 


III A834 


2.14el6 


III] A1664 


3.13el0 


IV A789 


1.17el6 


IV] A1401 


1.12ell 


VI A1034 


5.53el5 


Si III] A1892 


3.12ell 



Table 3 



a For T e = 10 4 K. 



Table 2: Analytic Line Flux Ratios 



Ratio 


B 


C 


Q 


-F 1 /F 2 a - 
= i q = n^i 


N 11} /C II] 


0.33 


-0. 


.53 


0. 


.045 


0.14 


N III]/C III] 


0.36 


-0. 


.69 


0. 


.060 


0.11 


n ni]/o in] 


2.19 


+0 


.42 





.36 


0.22 


N III] /Si III] 


0.067 


-0. 


.62 


0. 


,12 


0.007 


C III] /Si III] 


0.19 


+0 


.07 


1. 


.98 


0.065 


N IV] /O III] 


8.82 


+1 


.04 


0. 


.56 


0.90 


N IV] /O IV] 


4.58 


+0 


.59 


0. 


.85 


0.39 


N IV] /C IV 


0.28 


-0. 


.39 





,056 


1.4e-6 


N III/C III 


0.30 


+0. 


.21 


0. 


.11 


0.055 


N III/O III 


1.80 


+2. 


.73 


1 


,40 


0.53 


N V/C IV 


0.95 


-2. 


.32 





,087 


0.15 


N V/O VI 


1.09 


+2. 


.32 


0. 


.39 


0.25 


N V/(C VI+O VI) 








0. 


,071 


0.094 



a The flux ratios assume f{X\)/f(X2) — 1 and solar 
abundances, plus T 4 = 0.8 for N II}/ C II], T 4 = 2.2 for 
the N V ratios, and T 4 = 1.5 for all others. 



Model Abundances" 



z h 


N/O 


N 


He 


0.2 


0.2 


0.04 


0.94 


0.5 


0.5 


0.25 


0.96 


1.0 


1.0 


1.0 


1.00 


2.0 


2.0 


4.0 


1.07 


5.0 


5.0 


25.0 


1.29 


10.0 


10.0 


100.0 


1.66 



"All quantities are rela- 
tive to solar values. 
Z is the metals scale 
factor, £ (see Appendix). 
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Depth (cm) 



